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Objectives 


The  main  goal  of  this  program  has  been  to  support  the  immediate  needs  of  the  Air 
Force  as  well  as  being  relevant  to  long  range  Air  Force  defense-after-next  plans.  The 
program  includes  basic  research  that  is  expected  to  serve  a  role  in  the  long-range  goals 
of  the  Air  Force.  The  focus  of  this  research  has  been  to  support  the  scientists  at  Air 
Force  Research  Laboratory  at  Wright-Patterson  Air  Force  Base  involved  in  advanced 
technological  materials  processing.  The  foremost  goal  was  to  model  the  pulsed-laser 
ablation  plume  dynamics.  This  is  the  necessary  first  step  towards  a  process  modeling 
and  control  of  the  Pulsed-Laser-Deposition  (PLD)  process  of  materials  such  as  the  the 
high-temperature  superconductor  Yttriumi  Barium2  Copper  Oxygen  (YBCO).  The 
AFRL  has  a  multi-directorate  team  led  by  J.  Maguire  and  R.  Biggers  which  is  focussed 
on  making  large  specimens  of  YBCO.  However,  besides  YBCO,  there  are  other  related 
compounds:  CuO,  Ce02,  Zr02,  Y2O3,  and  YSZ.  Some  of  these  compounds  are  needed 
in  the  processing  of  YBCO,  while  others  are  simplified  material  models  to  aid  in  un¬ 
derstanding  YBCO.  Films  made  by  PLD  include  superconductors,  ferroelectric  ceramic 
films,  transparent  conducting  films,  lubricant  films,  and  semi-conductors.  Film  quality 
can  be  extremely  high;  however,  consistency  in  results  still  falls  far  short  of  what  is 
needed  before  turning  the  task  over  to  a  manufacturing  team.  The  main  factor  ham¬ 
pering  this  transfer  has  been  a  lack  of  model-based  control.  A  long-term  goal  is  the 
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a  matrix  which  may  or  may  not  itself  be  magnetic.  Most  of  the  quantitative  work  on 
the  magnetocaloric  effect  has  not  been  very  accurate.  This  is  because  the  interactions 
of  the  spin  thermodynamic  system  with  the  underlying  elastic  vibrations  has  not  been 
taken  into  account.  The  full  interaction  is  very  difficult,  since,  even  under  the  assump¬ 
tion  of  linear  elasticity  and  linearized  magnetism,  where  the  excitations  are  phonons 
and  magnons,  the  magneto-elastic  interaction  is  inescapably  non-linear.  An  objective 
of  this  research  program  has  been  to  study  models,  where  the  magnon  and  phonon 
branches  of  the  spectrum  combine  to  give  a  new  type  of  excitation  that  is  used  for  the 
quantum  statistics.  What  is  very  interesting  is  that  these  quasi-phonons  can  have  a 
spectral  branch  that  disappears  near  a  phase  transition. 


Accomplishments /Findings 


Concerning  the  ablation  plume  phenomenon,  a  computational  model  for  the  inter¬ 
action  of  laser  energy  with  solid  matter  was  developed,  as  well  as  analytical  and  com¬ 
putational  models  for  the  plume  hydrodynamics,  including  the  chemical  reactions.  A 
first  task  was  to  understand  which  variables  are  the  most  important.  There  are  approxi¬ 
mately  20  input  variables  that  can  be  varied  during  the  process.  Some  of  these  are  fixed 
parameters,  other  depend  on  time,  while  others,  such  as  electromagnetics,  are  (time- 
varying)  fields.  Contrasted  with  the  many  input  variables  are  the  desired  outcomes  for 
the  film.  Naturally,  there  is  not  time  for  experiments  that  vary  all  of  the  inputs  and  to 
measure  the  effect  on  these  outcomes  as  the  dimensional  complexity  is  overwhelming. 
We  needed  to  reduce  the  experimental  effort,  via  straightforward  polynomial  regression 
experimental  design  and  adaptive  learning  data  mining,  response  surface  methods,  The 
interaction  of  laser  energy  with  a  solid  material  takes  many  different  forms  depending 
on  the  power  input  and  the  material  properties.  We  have  modeled  this  using  a  phase 
change  approach  to  the  thermal  development.  That  is,  thermal  conduction  is  the  main 
method  of  energy  transport.  The  justification  is  that  the  energy  transferred  is  spread 
out  over  a  vibrationally  long  time  period.  The  important  thing  is  to  know  the  compo¬ 
sition  (in  multi-component  phase  space)  of  the  plasma  cloud  as  it  begins  its  expansion. 
There  is  an  important  role  played  by  the  hydrodynamics  (shock  heating,  mixing)  in 
the  stoichiometry  and  equilibrium  flux  onto  the  substrate.  A  focus  of  the  present  effort 
was  on  figuring  out  exactly  what  the  time-of-fiight  measurement  is  telling  us  about  the 
physics  in  the  plume  and  energy  flow  thereinto,  since  it  is  the  major  output,  and  is 
currently  used  as  a  manual  control.  There  are  interesting  features  of  the  time-of-flight 
signature,  such  as  double  peaking,  that  are  presently  the  object  of  intense  debate.  We 
have  identified  several  sources  for  these  phenomena,  including  reflected  shocks,  shock- 


complete  modeling  and  theory  necessary  for  an  eventual  automatic  control  of  this  pro¬ 
cess.  Towards  that  long  term  goal,  a  necessary  first  step  and  an  immediate  goal  of 
the  research  program  has  been  the  physical  and  mathematical  modeling  of  the  react¬ 
ing  plume  itself.  There  are  three  main  physical  subprocesses  occuring  during  the  PLD 
process.  These  are  the  1)  Ablation  Process,  where  photons  are  absorbed  by  the  target, 
resulting  in  heating  and  chemical  disassembly,  and  release  of  vapor,  2)  Plume  Dynamics, 
where  a  plasma  cloud  expands  and  suffers  self-collisions  and  collisions  with  the  ambient 
gas,  3)  Deposition  and  Growth.  Since  the  output  of  each  sub-process  supplies  input 
for  the  next  process,  there  has  been  some  expectation  that  the  modeling  should  follow 
this  sequence;  however,  we  have  found  that,  although  there  is  some  separation  of  time 
scales  from  one  sub-process  to  the  next,  there  is  also  considerable  feedback  and  mutual 
exchange  of  energy.  The  science  involved  is  at  the  extremes  of  energy  and  stability,  gov¬ 
erned  by  highly  nonlinear  processes,  mainly  due  to  the  interaction  of  disparate  physical 
scales  in  the  same  problem.  The  plume  is  very  dynamic  (supersonic,  high-temperature, 
shocks,  etc.)  especially  in  the  presence  of  a  reacting  background  gas.  (This  is  the 
case  for  deposition  of  YBCO,  which  requires  an  oxygen  ambient).  R.  Diggers’  lab  has 
been  running  depositions  for  several  years,  and  have  accumulated  much  data  concern¬ 
ing  the  plume:  spectroscopic,  time-of-flight,  spatial  distribution,  etc.  At  the  same  time, 
post-processing  measurements  have  been  made  of  the  resulting  film  (material  and  elec¬ 
tromagnetic).  A  typical  data  file  contains  over  80  different  entries  for  each  sample.  An 
objective  was  to  study  this  data  and  make  recommendations  whenever  necessary  and 
look  for  ways  to  use  the  information  in  model  building. 

If  the  pulsed  laser  deposition  process  is  ever  to  manufacture  films  that  are  needed 
by  the  Air  Force  (besides  many  other  material  systems  of  interest  to  AF  and  studied 
at  AFRL,  YBCO  coated-conductors  are  envisioned  for  use  in  generators,  especially  in 
pulsed-power  weapons  and  on  aircraft,  as  well  as  in  more  advanced  propulsion  systems), 
the  process  will  need  to  automatized.  Automatic  feedback  control  is  a  necessary  part  of 
this  process. By-products  of  the  same  research  program  will  be  better  material  quality, 
more  efficient  use  of  facilities,  and  better  system  throughput,  given  that  science  and 
mathematical  based  processing  control  is  easily  adapted  to  the  evolution  of  materials 
science  and  needs  of  the  Air  Force. 

Another  area  of  research  that  is  of  relevance  to  the  Air  Force  and  to  the  aerospace 
industry  is  modeling  of  magnetoelasticity,  and  more  particularly,  the  magnetocaloric 
effect  in  magnetic  and  superconducting  materials,  taking  into  account  magneto-elastic 
interactions.  The  magnetocaloric  effect  has  been  used  at  low  temperatures  in  adiabatic 
demagnetization,  but  recently  has  been  in  the  news  as  a  possible  method  for  refriger¬ 
ation.  A  prototype  of  a  refrigerator  using  the  magnetocaloric  effect  has  been  built  at 
the  Astronautics  Corporation.  The  materials  used  in  these  magnets  are  superparamag- 
netic  nanoclusters:  that  is,  they  consist  of  magnetic  nano-sized  particles  imbedded  in 


induced  vorticity,  charge  separation  and  chemical  species  separation.  There  are  reports 
from  both  experiment  and  simulation  that  show  exagerated  forward-direction  to  the 
plume.  We  have  an  analytic  model  to  explain  this.  We  have  also  been  using  simplified 
material  systems  (Cu,  CuO,  CeO)  as  tests  of  the  basic  model.  Based  on  numerical  ab¬ 
lation  models,  the  number  density  of  an  ablated  plume  ranges  from  1016  cm-3  to  above 
1020  cm-3,  depending  on  the  energy  input.  With  the  energies  used  in  film  deposition 
in  a  background  gas,  one  sees  the  latter  range.  Total  particle  energies  in  the  nascent 
plume  are  typically  in  the  range  4-17  eV.  The  translational  energy  is  thermalized  af¬ 
ter  only  a  few  collisions  per  particle;  we  estimate  a  Knudsen  number  on  the  order  of 
10"4  —  10-2  in  the  flow  domain  behind  the  shock  front  -  a  mean  free  path  is  estimated 
at  <  0.1  mm.  All  the  same,  the  particles  evaporating  or  being  ejected  from  the  target 
will  have  a  non-equilibrium  distribution  that  will  relax  to  a  Maxwell-Boltzmann  dis¬ 
tribution  after  several  (3-10)  collisions.  For  ablation,  or  beam  evaporation  from  a  hot 
chamber,  this  takes  place  in  a  one-dimensional  region  nominally  <  1  mm  thick.  For 
this  reason,  we  treat  this  region  here  as  an  infinitesimal  front  much  like  a  stationary 
shock.  For  given  material  parameters  and  laser  beam  energy,  this  and  various  other 
simulations  and  calculations  supplies  us  with  the  values  of  hydrodynamic  variables  at 
the  “outer”  edge  of  the  Knudsen  layer,  which  is  the  inlet  region  for  large-scale  numerical 
computations.  In  very  low  pressure  conditions  the  plume  expansion  can  be  described 
by  a  adiabatic  expansion  into  vacuum.  For  expansion  into  a  higher  pressure  ambient, 
there  are  two  competing  models  for  one  dimensional  plume  expansion.  The  blast  wave 
model  describes  the  self-similar  radial  expansion  when  energy  is  deposited  at  a  point  at 
t  —  0.  The  normal  shock  model  describes  one-dimensional  motion  peaked  in  a  partic¬ 
ular  direction.  Following  closely  on  experiments,  special  solutions  can  be  constructed 
each  valid  on  scales  widely  separated  from  others.  Computations  also  show  this  dif¬ 
ference.  However,  in  some  situations  the  plume  dynamics  cannot  be  described  by  a 
one-dimensional  model.  The  plume  dynamics  in  the  presence  of  a  temperature  gradient 
is  one  such  example.  Particularly  in  situations  when  the  plume  is  seen  to  be  extremely 
forwarded-directed.  In  experiments  where  substrate  heat  was  isolated  as  a  control,  the 
plume  showed  an  increase  in  size  and  emission  intensity  that  was  similar  to  the  effect 
of  increasing  laser  beam  energy.  However,  in  a  vacuum,  the  heater  produced  no  dis¬ 
cernible  change  to  the  plume.  Thus,  the  effect  was  due  to  molecular  kinetics  and  not  to 
radiative  heat  transfer.  Moreover,  this  has  been  reproduced  in  our  numerical  computa¬ 
tions,  and  has  been  explained  via  an  analytical  model.  Two  dimensional  calculations  of 
the  complete  multi-component  reacting  continuum  fluid  were  carried  out  under  a  wide 
range  of  operating  variables,  including  oxidation  reactions.  The  main  features  of  the 
plume  shape  in  the  presence  of  substrate  heating,  however,  can  be  demonstrated  with  a 
one-component  fluid.  The  difference  between  the  computed  plumes  without  substrate 
heating  differ  dramatically  from  the  computed  plumes  in  the  presence  of  a  temperature 
gradient,  and  corresponds  to  the  difference  that  one  sees  in  actuality.  The  development 
of  vorticity  is  clearly  indicated,  as  well  as  the  one-dimensional  jet-like  flow.  A  compre- 


hensive  collection  of  simulation  results,  including  animations,  can  be  found  on  the  web 
site  of  the  PI  listed  above. 

The  excimer-type  laser  pulse  lasts  between  tens  of  femtoseconds  up  to  hundreds  of 
nanoseconds.  During  this  time,  which  can  be  long  relative  to  the  times  characteristic  of 
the  basic  interaction,  the  target  material  can  become  multi-phased  as  the  solid  melts  and 
is  ablated,  either  evaporatively  or  explosively.  Furthermore,  there  is  sometimes  charge 
separation  as  a  plasma  cloud  forms.  Thus,  there  are  several  competing  pathways  for 
the  light  energy  to  be  converted  to  other  forms.  On  the  shortest  time  scales,  that  of  the 
optical  excitation  of  electrons  in  a  metal  or  of  excitation  of  quasi-particles  in  a  dielectric, 
one  must  take  into  account  the  disparate  relaxation  times  of  the  electron  system  and 
the  lattice-vibrational  system  of  the  solid.  Thus,  we  consider  two  temperatures,  Te, 
the  temperature  of  the  electrons,  and  Ti,  the  temperature  of  the  ions  in  the  lattice. 
The  interaction  between  these  thermodynamic  systems  can  be  modeled  via  diffusion 
equations  with  a  weak  coupling  between  them: 

c‘f  =  |wr-)f)+B(T--T‘>’ 

where  Q,  the  source  term,  is  the  laser  power  absorbed  (in  this  example  by  electrons, 
bound  or  in  conducting  band),  and  2  is  the  coordinate  normal  to  the  surface.  It  is 
reasonable  to  assume  that  the  process  is  one-dimensional  on  the  fast  time  scales  where 
this  model  is  valid.  At  later  stages  (on  longer  time-scales),  the  situation  may  be  higher 
dimensional,  but  by  that  point  the  two  thermodynamic  systems  can  be  considered  to 
have  achieved  equilibrium  with  one  another.  Higher  dimensional  models  are  only  really 
necessary  for  long  pulsed  systems  or  continuous  (CW)  laser  work.  The  correct  scaling 
depends  on  the  the  power  absorption;  for  high  power  absorption  (~GW  cm-2),  the 
two  temperature  model  is  necessary  in  order  to  model  the  ablation,  but  for  moderate 
intensities,  where  power  is  delivered  over  tens  of  nanoseconds,  a  single-temperature 
multi-phase  thermal  model  may  be  sufficient.  The  material  properties,  like  the  optical 
properties,  depend  generally  on  the  temperature  and  may  be  discontinuous  at  a  phase 
boundary,  so  that  we  have  the  pair  of  equations, 

c*<T>4r-  =  +<3-  (1) 

“<r>4r‘  -  lk'm §+«-•  <2> 

where  the  subscripts  s  and  l  refer  to  the  solid  and  liquid  regions,  respectively.  These 
equations  are  supplemented  by  the  boundary  condition, 


Ts  To,  z  — >  00, 


(3) 


where  T0  is  the  controlled  heat  of  the  workpiece  or  target,  as  well  as  conditions  at  the 
phase  interfaces,  z  =  asi(t),  and  2  =  aiv(t).  At  the  solid/liquid  interface,  conservation 
of  energy  requires 

,  dTs  dT,  _  d 
ks  dz  kl  dz  ~ PsLm dta$u 

where  Lm  is  the  latent  heat  of  fusion.  Superheating  has  been  shown  to  have  thermal 
effects  on  the  order  of  several  percent  for  some  laser  processing  applications,  however  for 
the  nanosecond  timescales  which  which  are  our  main  focus  it  is  considered  negligible. 
Most  important  for  high  energy  applications  (such  as  ablation)  is  the  possiblity  of  non- 
equlibrium  processes,  such  as  phase  explosion.  Such  modeling  is  part  of  the  current 
proposal.  It  is  a  straightforward  task  to  verify  the  consistency  of  this  assumption. 
Energy  conservation  at  the  liquid/vapor  interface  requires 

,zm_  r  £ 

kldz  -  plL* dtaiv’ 


and  the  kinetics  of  evaporation  mean  that 


dh  =  ~7,aiv  =  v*e  Ah/kTlv, 
at 

where  A h  is  enthalpy  of  vaporization  (per  molecule),  vt  is  approximately  equal  to  the 
speed  of  sound  in  the  material,  and  is  the  temperature  at  the  interface.  A  result 
of  these  relations  is  the  fact  that  Tiv  can  be  much  larger  than  Tv,  the  normal  boiling 
temperature.  This  temperature  also  depends  on  the  vapor  saturation  pressure.  The 
laser-induced  heat  flux  is  given  by 

9  =  (1  -  R(T))e-f°°<T(*'»d/I(x,y,t), 


where  R  is  the  reflectance,  which  depends  on  the  photon  frequency  as  well  as  on  the 
temperature  and  phase  of  the  material,  a  is  the  linear  absorbtion  coefficient  (which 
also  depends  on  photon  frequency  and  temperature  and  phase  of  the  material),  I  is 
the  intensity  of  light  radiation  incident  on  the  surface,  and  z  is  the  distance  into  the 
condensed  material,  normal  to  the  surface.  The  heating  rate  density  source  term  is 
thus,  for  constant  absorption  coefficient, 


Q  =  (1  -  R)ae  azI(x,y,t). 


The  laser  spot  is  not  symmetric  and  generally  a  non-trivial  function  of  the  tangential 
variables,  x  and  y.  Although  the  spatial  profile  of  the  excimer-type  laser  illumination 
spot  is  usually  taken  to  be  roughly  rectangular,  it  varies  with  time  and,  particularly, 


between  shots,  in  unpredictable  ways.  The  fact  is,  is  it  very  difficult  to  measure  irradi- 
ance  at  all,  let  alone  in  real  time.  The  laser  pulse  intensity  will  be  assumed  to  have  the 
characteristic  excimer  form: 

I(t)  =  /0(—-)Me~</Tp, 
rv 

rather  than  the  Gaussian  pulse  that  is  sometimes  used  (and  characteristic  of  the  older 
generation  of  lasers).  The  relation  that  rp,  which  we  will  call  “pulse  length,”  has  to  the 
full-width  half-maximum  (FWHM)  of  pulselength  is 

tp/tf  «  1.2. 


There  are  several  approaches  to  the  scaling  of  the  problem.  If  we  take  as  charac¬ 
teristic  time  scale  for  the  problem,  the  pulselength  tp,  and  from  this  form  the  diffusion 
length,  Id  —  ^ Dtp ,  which  we  take  to  be  the  characteristic  length  scale.  The  diffusivity, 
D,  is  defined  to  be 


ks(To) 

cs(To)ps 


(4) 


If  material  parameters  are  assumed  to  be  independent  of  temperature  within  each 
phase,  then  the  equations  become: 

(5) 

where 

/,((,r)  =  e'(oW'“,')r'‘e‘r.  (6) 

The  dimensionless  parameter, 

_  q(l  -  R(Tm))I0Tp 
cp(Tv  -  T0)  ’ 

appears  multiplying  the  source  (6).  From  the  form  of  the  equation  (5),  it  follows  that 
the  dimensionless  number,  V/alp  is  a  measure  of  the  ratio  of  absorbed  light  energy 
density  to  energy  density  needed  to  reach  the  vaporization  temperature.  Numerical 
computations  bear  out  the  intuition  that  T  is  indicative  of  maximum  temperature  rise, 
while  the  dimensionless  group  T  joio  gives  a  lowest-order  prediction  of  the  overall  phase 
change. 

The  boundary  condition  (3)  becomes 


6  — t  0,  C,  — y  oo. 


A  purely  thermal  condition  at  the  melt  interface  is 

d6s  dO  i  da 

where  ft  =  ks/ki  is  the  ratio  of  the  conductivities  at  the  melt  temperature,  and  rjm 
is  a  dimensionless  constant  characterizing  the  relative  size  of  the  enthalpy  of  fusion 
compared  to  the  rest  of  the  energy  required  to  reach  the  boiling  temperature  (this 
number,  for  constant  parameter  values,  is  the  reciprocal  of  the  Stefan  number).  Values 
of  fd  are  typically  ~  2.25  —  2.5  for  metals,  «  1  for  ceramics. 

The  enthalpy  method  has  been  used  for  computations  because  it  can  easily  incorpo¬ 
rate  the  nonlinearities  due  to  the  temperature  dependencies  of  the  material  parameters, 
and,  since  we  are  mainly  concerned  with  the  approximate  size  and  evolution  of  the  melt 
region,  we  don’t  have  to  worry  about  the  usual  difficulties  with  front-tracking.  The 
enthalpy,  or  heat  content  per  unit  volume,  is  given  by 

H=  \T  pc(r)  +  pLm6(T'-Tm)dT', 

J  0 

where  the  specific  heat  follows  an  empirical  law, 

c(T)  =  a  +  6T  +  cT2  +  dT 3  +  eT~2, 

where  the  coefficients  are  tabulated  in  various  references  or  else  matched  to  data  tables 
or  graphs  from  other  sources.  The  temperature  dependency  of  the  thermal  conductivity 
is  taken  care  of  using  the  Kirchoff  transformation, 


The  nonlinear  function,  k(T),  is  much  harder  to  match  to  a  general  formula  in  integral 
powers  of  the  temperature.  A  weak  formulation  of  the  model  is  then  considered  using 
the  variable  H  and  W  by  using  an  integral  operator  form  of  the  heat  conservation  law. 
This  results  in  a  toughening  of  the  the  usual  stability  condition  for  the  explicit  method, 
and  in  practice,  the  time-step  is  necessarily  much  smaller  than  order  of  square  of  the 
space  interval. 

Figure  2  portrays  the  results  of  a  calculation  using  material  parameters  typical  of  a 
“black”  oxide,  with  optical  properties  between  those  of  a  semi-transparent  oxides  and 
metals.  A  small  absorption  coefficient  and  low  thermal  conductivity  combine  to  ensure 
that  the  material  “sees”  a  practically  uniform  body  heating.  Physically  this  should 
correspond  to  the  establishment  of  a  “mushy”  region.  As  is  seen  in  Figure  3,  the  rise 
to  maximum  temperature  (just  above  the  melting  temperature)  is  fast.  The  computed 


“mushy”  region  is  clearly  visible  and  lasts  several  pulselengths;  however,  when  kinetics 
of  nucleation  and  interface  formation  are  taken  into  account,  along  with  the  minimal 
temperature  rise  over  the  melting  temperature,  one  expects  the  physical  “mushy”  region 
to  be  more  considerable.  The  time  scale  is  about  30  pulse  lengths  (close  to  1  /is). 

In  nanosecond  pulses  of  intensity  >  108  W  cm-2,  energy  is  transferred  to  the 
solid/liquid  at  such  high  heating  rate  that  there  is  a  critical  or  supercritical  transi¬ 
tion  through  a  metastable  phase  of  solid/liquid/gas.  In  this  regime,  the  material  is 
far  from  equilibrium,  as,  thermodynamically  it  can  be  thought  of  as  in  a  state  near 
the  spinodal  line  where  volume  fluctuations  are  high.  This  means  that  the  material 
undergoes  homogeneous  nucleation  of  vapor  regions  inside  the  condensed  phase  rather 
than  usual  boiling.  The  energy  may  be  released  in  a  shock  wave  of  the  multi-phase 
condensed  material.  The  amount  of  ablated  material  depends  linearly  on  the  fluence 
of  the  laser  beam,  as  long  as  the  intensity  is  not  too  high  or  too  low.  The  ablated 
material  is  at  extremely  high  temperatures  (~  104  —  105K)  and  thus  ionized;  the  degree 
of  ionization  varies  according  to  circumstances,  and  has  been  measured  as  being  much 
higher  than  predicted  by  equilibrium  considerations  (Saha  equation).  Besides  the  ther¬ 
mal  ionization  of  the  plume,  there  are  other  pathways,  including  photo-ionization  and 
multiphoton  absorption.  The  ionized  cloud  erupting  from  the  surface  can  absorb  energy 
from  the  laser  beam.  Thus  the  incipient  plume  is  further  heated  by  this  absorption.  It 
is  expected  that  the  dominant  means  of  photon  absorption  in  this  gas  is  by  three-body 
interaction  involving  two  electrons  and  a  neutral  atom;  the  degree  of  ionization  of  the 
plume  has  never  been  conclusively  established  by  experiment  and  is  presumably  very 
sensitive  to  control  parameters.  If  the  ablation  is  done  into  vacuum,  the  expansion  is 
rapid  enough  that  the  bulk  of  the  plume  is  transparent  to  the  laser  beam,  and  the  ab¬ 
sorption  region  is  at  the  gas/condensed  interface.  This  layer  is  very  small  relative  to  the 
size  of  the  expanding  plume  and  can  be  considered  to  be  held  at  a  single  temperature, 
since  the  relaxation  here  is  much  faster  than  the  accelerating  plume.  Expansion  into 
a  vacuum  proceeds  according  to  free  flight  where  the  thermalization  of  particles  never 
takes  place.  For  near- vacuum  depositions,  the  flight  of  the  plume  can  be  expected  to 
be  a  Knudsen  flow  where  collisions  slowly  establish  local  thermodynamic  equilibrium. 
Film  depositions  are  usually  done  in  an  oxidant  ambient  of  about  50  —  300  mTorr. 
In  this  atmosphere,  the  ablation  material  is  an  accelerating  jet  into  the  ambient  and 
a  compressive  shock  is  formed.  This  forms  the  supersonic  plume.  A  one-dimensional 
model  has  been  constructed  using  the  assumption  of  equlibrium  phase  kinetics  and  the 
use  of  the  Knudsen-Layer  as  a  hydrodynamic  discontinuity.  For  this  model,  the  result¬ 
ing  hydrodynamics  depends  on  the  single  parameter  of  the  surface  temperature,  which, 
however,  depends  on  control  parameters  in  a  complicated  way. 

Regarding  the  hydrodynamics  and  chemistry  of  the  plume,  we  developed  both  ana¬ 
lytical  and  computational  models.  We  were  particularly  interested  in  resolving  questions 


that  arose  during  analysis  of  experimental  data,  particularly  time-of-flight  measure¬ 
ments  from  emission  spectroscopy.  In  the  plume  itself  the  general  conservation  laws  are 
valid.  That  is,  for  general  regions,  V ,  and  times,  a  system  of  integral  equations, 

/  G(r,t2)-G(r,tl)dr=  [  F  •  dS, 

Jv  JdV 

holds,  where  G  is  a  vector  of  conserved  quantities  and  F  is  a  tensor  of  fluxes.  In  order 
for  this  to  reduce  to  a  differential  model,  we  need  to  identify  some  small  size  V  so 
that  thermal  equilibrium  (among  any  components  lumped  together)  is  instantaneously 
established  on  this  scale  (as  far  as  the  model  is  concerned).  Here,  we  identify  such  a 
dimension  with  the  mean  free  path  based  on  densities  encountered  in  the  deposition 
chamber.  The  pulsed  laser  ablation  plume  is  made  up  of  M  components,  of  varying 
mass.  Let  n,-  denote  the  (particle)  concentration  of  the  ith  species,  and  mi  the  (particle) 
mass  of  the  ith  species,  then  define  the  peculiar  velocities  in  the  following  way:  each 
particle  of  the  ith  species  has  a  velocity 

Vi  =  v  +  Vi 

where  v  is  the  bulk  flow  velocity  and  V,-  is  the  peculiar  velocity  of  the  ith  species.  Note 
that  in  the  particle  phase  space,  each  particle  takes  on  the  local  momentum: 

miVi  =  p. 

After  Boltzmann  averaging  and  mass-averaging  over  the  species,  we  get  a  local  flow 
velocity  defined  by 

pv  =<  X>v,  >p=  ^pi  <  Vi  >p, 

where  <  >p  denotes  the  Boltzmann  averaging  over  momentum  space.  The  peculiar 
velocity  has  the  average, 

<  ^  >p=  yf, 

which  we  call  the  molecular  diffusion  velocity,  so  that 

Vi  =  v  +  Vf  +  V-, 

where  V(-  is  the  “chaotic  motion”  of  the  ith  species.  This  allows  us  to  define  a  temper¬ 
ature  for  the  ith  species, 

T-  =  — 

2c„i 

By  averaging  over  scales  smaller  than  the  characteristic  length  scale  identified  above,  we 
can  derive  the  conservation  laws  valid  on  larger  scales.  The  species  mass  conservation 
takes  the  form 


<  iv;i2  > 


|ni  +  V-(ni(v  +  Vf))  =  5i,  1  =  1,...,M, 


where  p,-  is  the  species  production  or  decomposition  due  to  chemical  reactions.  This 
equation  could  also  be  written  in  terms  of  the  mass  densities,  p,  =  n.m,-,  or  the  mass- 
fractions,  pi  I p.  Multiplying  by  m,  and  summing  over  i,  and  using 

^mtn;  <  Vf  >=  0, 

we  get  the  bulk  mass  conservation  equation, 

JtP  +  V  •  (pv)  =  0. 


The  equation  of  motion  for  the  mixture, 

d 

p-^y  +  pv  •  Vv  +  V  ■  V  =  0, 

is  written  in  terms  of  a  total  pressure  tensor.  If  we  assume  that  <  V;  >2<C<  Vj2  >, 
which  will  be  valid  except  in  “collisionless”  flow,  then  we  can  derive  equations  for  the 
.  temperatures.  The  total  pressure  tensor  is  written 

?  =  £«, 

where  the  pressure  for  each  component  is  that  of  an  ideal  gas, 

Rt 

Pi  =  Pi — Ti. 
rrii 

The  overall  mixture  will  then  satisfy  an  ideal  gas  law, 


p  =  nRT, 


if  we  define  a  (mixture)  temperature  as 


r  =  -£»*, 


n 


where  n  =  £n,-.  (This  was  the  whole  point  of  defining  things  this  way.  We  will  deal 
with  an  inviscid  fluid,  so  there  is  no  ambiguity  about  the  meaning  of  pressure.) 


The  specific  enthalpy, 

h  =  J  CpdT  +  p-1  ^p,- A/ifi, 

is  conserved.  Here,  A jH{  is  the  heat  of  formation  of  the  ith  species,  and  T  is  the 
total  mixture  temperature  as  defined  above.  In  this  way,  the  energy  source  (sink)  of 


an  exothermic  (endothermic)  reaction  is  taken  into  account.  At  the  rapid  flow  acceler¬ 
ations  in  the  plume,  the  entropy  variation  is  small  and  therefore  the  reaction  term  is 
immediately  converted  into  a  change  in  temperature. 

It  is  easy  to  estimate  characteristic  scales  and  relative  importance  of  transport  mech¬ 
anisms,  based  on  a  one-dimensional  model  of  the  preparation  of  the  blow  up  to  the 
nether  edge  of  the  Knudsen  Layer.  A  mean  free  time  between  collisions  is  estimated  as 
10  —  20  ns,  while  the  characteristic  time  for  total  (mean)  momentum  to  be  transferred 
across  features  of  the  plume  by  collisions  is  about  100  times  longer.  This  is  the  main 
reason  that  we  consider  an  inviscid  flow.  However,  because  of  the  large  temperature 
gradients  that  are  set  up  by  the  shock  heating  in  the  supersonic  plume,  and  the  in¬ 
teraction  with  chemical  reactions,  we  must  allow  for  thermal  conduction  and  species 
diffusion. 

In  short,  the  inviscid  hydrodynamic  equations  with  heat  and  mass  diffusion  and 
heat  production  by  reactions,  comprise  our  basic  model  according  to  the  approximation 
above.  We  carry  out  a  preliminary  one-dimensional  inviscid  calculation  and  then  a 
two-dimensional  transport  calculation  with  the  chemistry  accounted  for  in  a  simplified 
scheme  which  is  explained  in  what  follows,  as  we  focus  on  the  chemistry  of  the  plume 
(Figure  2).  Metal  vapor,  atomic  and  molecular  oxygen  (and/or  other  oxidizers)  and  any 
relevant  metal-oxide  are  the  components  of  the  plume.  We  will  make  the  case  that  the 
primary  reaction  in  our  scheme  is  the  termolecular  oxidation  reaction.  In  the  nascent 
plume,  one  expects  to  see  bimolecular  exhange  reactions,  governed  by  arbitrary  high 
kinetic  energy  collisions,  such  as 

M  +  02  — *  MO  +  0 

These  reactions  have  high  activation  energies  and  should  constitute  the  main  scheme 
in  the  Knudsen  layer  and  at  its  edge.  After  a  period  that  is  on  the  order  of  10  —  100 
ns,  the  metallic  vapor  has  been  excited  and  starts  recombination  and  de-excitation  due 
to  collisions  and  photon  release.  The  characteristic  times  of  the  bimolecular  reactions 
show  the  scaling 

Tbimol  OC  pke~Ea/kT. 

These  “endothermic”  reactions  establish  strong  localized  temperature  gradients.  On 
the  time  scales  of  the  fully  formed  plumes  simulated  here,  termolecular  reactions,  such 
as  the  oxidation 

2M  +  02  — >  2MO, 

or  the  recombination  reaction 


will  predominate.  These  are  heat  producing  reactions  that  are  self  sustaining  on  rel¬ 
atively  long  time  scales,  even  as  the  plume  is  cooling  due  to  expansion.  For  these 
reactions,  the  characteristic  time  scale  is 

Ttermol  PM") 

proportional  to  the  density  of  metal  atoms,  and  we  expect  that  this  will  lead  to  equilib¬ 
rium  densities.  In  this  way,  we  are  led,  in  short  order,  through  an  initial  oxide  formation 
via  the  endothermic  bimolecular  reactions  in  the  nascent  plume,  cooling  by  expansion  as 
the  plume  loses  its  one-dimensional  character,  balanced  by  the  heat-producing  reactions 
which  are  sustained  through  until  equilibrium  is  achieved. 

Figure  4  shows  the  results  of  a  two-dimensional  simulation  of  the  plume  equations 
carried  out  by  use  of  a  finite-volume  method.  The  top  row  shows  the  temperature 
distributions  at  2  different  times,  (a)  0.5/rs,  (b)  1.0/is,  (c)  1.5ps,  for  a  case  of  medium- 
high  energy  input.  The  second  row  shows  the  density  and  pressure  distributions  for  the 
same  case.  The  third  row  shows  the  density,  temperature,  and  pressure  distributions 
for  a  case  for  smaller  energy  input  at  a  time  of  2/xs.  The  approximate  size  of  each 
picture  is  a  centimeter.  Thus,  the  lower  energy  jet  has  a  less-pronounced  shock,  and 
more  spherically  symmetric  temperature  distribution. 

The  difference  between  the  computed  plumes  without  substrate  heating  (Figures  4 
and  5)  differ  dramatically  from  the  computed  plumes  in  the  presence  of  a  temperature 
gradient  (Figures  6,  7,  8),  conforming  to  the  difference  that  one  observes  in  the  actual 
plume. 

The  pictures  of  a  plume,  either  in  vacuum,  or  in  ambient  without  substrate  heating 
show  the  usual  rounded  cosine  shape,  but  when  substrate  heating  is  present  strong 
vorticity  is  generated  at  the  ‘wings’  of  the  plume,  leading  to  recirculating  flow,  greater 
mixing,  and  thermalization  in  those  regions,  but  allowing  a  more  focussed  quasi-one- 
dimensional  jet  to  be  produced  along  the  flow  axis  as  the  plume  nears  the  substrate 
heater.  One  of  the  prominent  features  observed  in  the  simulations  are  relatively  large 
scale  vortices  at  the  oblique  wings  of  the  plume.  Recirculating  flow  and  strong  mixing 
occur  there  giving  rise  to  collisional  processes  that  result  in  a  strong  emission  sometimes 
showing  a  double  peak  in  the  density  time-of-flight.  The  sense  of  the  vortices  is  what  we 
will  refer  to  as  expected  orientation.  A  comprehensive  collection  of  simulation  results, 
including  animations,  can  be  found  on  the  author’s  web  site. 

The  expression  for  the  vorticity,  u>  =  curlu,  is,  in  these  coordinates, 

_  ,dw  1  dijxjv).  [  t  1  d{hu)  dw  t  t  1  d(h2v)  1  d{hxu ) 
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Figure  4:  Computed  shock  wave  of  an  ablation  plume.  Temperature  contours  are  shown. 


Figure  5:  Computed  shock  wave  of  an  ablation  plume.  As  in  Figure  4,  temperature 
contours  are  shown.  The  operating  conditions  are  different  than  in  that  case. 


Figure  6:  Temperature  distribution  (color)  and  velocity  vectors  at  three  successive 
snapshots  of  a  plume  in  an  applied  temperature  gradient,  taken  about  1  /is  apart.  It  is 
apparent  that  in  this  case  a  strong  recirculating  flow  is  developing. 


Figure  7:  Contours  of  the  temperature  distribution  in  the  devloping  plume,  for  a  set 
of  operating  conditions  similar  to  Figure  5  except  for  the  presence  of  a  temperature 
gradient. 


Figure  8:  Enthalpy  (color)  and  velocity  vectors  in  yet  another  case.  This  case  corre¬ 
sponds  to  a  greater  energy  deposition. 


|  Wia n  ]  =  n  X 


(\  -]' VtanK)  +  S  p](  — )(S(uian  -  V,on)  -  n)^  ,  (17) 
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where  we  have  only  included  terms  that  are  consistent  with  the  assumption  that  the 
shock  moves  into  quiescent  gas. 

It  is  the  first  term  in  the  above  expression  that  involves  the  thermal  gradient  in 
the  quiescent  gas.  The  other  term,  involving  the  jump  in  density,  is  made  up  of  two 
terms,  the  first  is  related  to  the  curvature  of  the  shock  front,  and  the  second  involves  the 
local  rotation  of  the  shock.  In  order  to  estimate  the  relative  importance  of  these  terms 
an  analytic  expression  for  the  thermal  field  in  the  quiescent  gas  will  be  used.  We  can 
approximate  the  heated  substrate  by  a  flat  disk,  radius  a,  held  at  a  fixed  temperature, 
Ts ,  and,  given  the  geometry,  sizes,  and  materials  involved,  the  temperature  at  infinity 
can  be  taken  to  be  room  temperature,  Tr.  Using  the  nondimensionalization,  0  =  , 

and 

P  _  1  _ _ 1 

prtI  1+0  1  +  0sarccot(^)’ 

where  6S  is  a  nondimensional  parameter  of  order  approximately  unity  when  the  substrate 
temperature  is  in  the  range 

Tr  <  Ts  <  ATr , 

which  is  essentially  the  experimentally  considered  range.  It  is  reasonable  to  assume 
that  the  density  behind  the  shock  is  constant.  From  the  functional  form  of  (18)  it 
is  seen  that  there  are  essentially  two  regimes:  for  distances  from  the  target  on  the 
order  of  a,  the  term  j  is  at  least  an  order  bigger  than  |  p] .  Let  us  examine  this 
expression  locally.  Assume  a  planar  shock,  oriented  parallel  to  the  i  axis  where  the  k 
axis  is  perpendicular  to  the  substrate  pointing  at  the  target,  with  V  =  — V(t)k,  and 
n  =  —  cos  0k  +  sin0j  —  sin^>k  +  cos^j.  (Here,  is  the  complement  of  0.)  This  0  is  a 
local  parameter,  but  can  be  associated  with  the  global  variable,  0 ,  for  a  shock  that  is 
approximately  spherical.  (By  choosing  this  form  for  the  tangent  plane,  we  are  assuming 
an  axi-symmetric  plume.)  In  the  same  way,  we  will  use  the  variable,  z,  which  can  be 
replaced  by  in  a  global  treatment.  The  mass  flux  is  ms  =  psin^V,  and 


Vtan(m)  =  -cost/;  sin  (19) 

Oz 

Since  f  ^  ]  is  negative,  the  vorticity  produced  is  of  the  expected  orientation.  Furthermore, 
in  the  early  regime,  before  the  plume  decelerates  and  while  that  factor  is  large,  the  term 
(19)  is  small  along  the  longitudinal  axis  of  the  plume  and  only  is  significant  off-axis. 
Now  let  us  examine  the  rotation  term,  using  the  same  local  plane.  Here,  n  =  fi  x  n. 
First,  consider  the  rotation  towards  or  away  from  the  substrate,  i.e.,  fl  =  wi,  with  to  <  0 


or  u  >  0,  respectively.  This  gives  rise  to  a  transverse  vorticity  production  term, 

m 

Furthermore,  this  term  is  opposite  signed  to  that  which  would  be  produced  due  to 
the  curvature.  An  arbitrary  rotation  can  produce  streamwise  vorticity  as  well;  let 
fi  =  wxi  +  u>yj  +  w2k,  then  the  corresponding  vorticity  production  term  is 

Ip]  (  — )  (— wxi  -f  (u)z  cos ip-uiy  sin ip)X) , 
m 

where  X  is  a  unit  vector  along  the  tangent  plane,  perpendicular  to  i.  Note  that  stream- 
wise  vorticity  production  can  be  zero  even  in  the  presence  of  rotation.  This  generation 
of  vorticity  is  especially  important  when  one  takes  into  account  the  interaction  of  re¬ 
laxation  and  reaction  mechanisms.  This  is  the  focus  of  ongoing  research. 

Concerning  the  theory  of  magnetocaloric  effect  in  magnetoelasticity,  we  can  derive 
interesting  results  under  the  condition  of  weak,  but  non-negligeable  interaction.  This 
weak  collective  excitation  scheme  was  then  carried  over  by  analogy  to  the  excitations  of 
the  flux-line  lattice  in  type-II  superconductors  (HTSC)  in  the  mixed  state.  It  may  seem 
at  first  that  this  is  a  completely  different  situation,  but,  as  long  as  the  same  statistics 
hold  (this  is  regime  of  linear  excitations),  the  situations  will  be  similar  when  the  energies 
behave  similarly. 

The  history  of  investigations  into  the  magnetocaloric  effect  has  shown  that  sharp 
changes  in  the  isentropic  curves,  are  associated  with  phase  changes.  These  investigations 
usually  rigorously  focus  on  the  entropy  of  the  magnetic  system  and  do  not  take  explicit 
account  of  the  vibrations  of  the  material  structure,  except  to  accord  it  the  status  of  heat 
sink.  In  a  magnetic  system,  for  example,  there  is  a  coupling  between  the  excitations 
of  the  magnetic  systems  (magnons)  and  the  vibrations  of  the  crystal  lattice  (phonons) 
and  the  condition  becomes 

^mag  Sphon  —  constant. 

If  the  crystal  lattice  is  simply  a  heat  sink  then  the  effect  of  the  phonons  is,  broadly  put, 
to  weaken  the  magnetocaloric  effect,  but  does  not  qualitatively  affect  the  calculation,  in 
particular,  the  sign  of  the  effect.  The  criteria  governing  the  validity  of  this  assumption 
can  be  made  explicit  by  reference  to  the  relevant  relaxation  times  of  the  various  modes 
of  dissipation. 


where 


^magnetic  =  flA/2  +  A(VM)2  _  .  z)2  _  M  •  H, 

where  M  is  the  (macroscopic)  magnetic  moment  density,  H  is  the  external  magnetic 
field,  a,  A  are  exhange  coefficients,  /?  is  the  anisotropy  coefficient  (the  easy  axis  of 
magnetization  is  the  z-direction).  The  crystal  is  assumed  to  be  isotropic  (cubic),  and 
the  energy  stored  in  the  elastic  deformation  is 

Hdastic  =  2Xull  +  ^u„k  +  uk,{)\ 

where  u  =  (ux,  %,  uz)  is  the  material  deformation  vector, 

duk 

Uk>i~  dXi' 

and  X,fx  are  the  isotropic  elastic  moduli.  The  interaction  energy  is  given  by 

jjmagnetoelastic  _ 


The  Hamiltonian  was  diagonalized  by  a  method  of  second  quantization,  leading  to  a 
8t/l-order  secular  equation  for  the  energies  of  the  collective  excitations  of  the  combined 
magnetoelastic  system.  Special  solutions  can  be  found  corresponding  to  a  (quasi-) 
magnon  branch  and  three  (quasi-)  phonon  branches.  In  other  words,  we  continue  the 
branches  that  start  at  g  =  0,  corresponding  to  spin  waves  and  to  sound  waves  propa¬ 
gating  along  symmetry  axes  having  two  degenerate  transverse  waves.  This  spectrum  is 
valid  under  the  assumption  of  weak  interaction 

(A  1/(3  «  1). 

Then,  the  results  of  the  calculation  are  that 
The  quasimagnon  spectrum  is 

fl2,  =  4  +  w{ki  +  k]). 


The  quasiphonon  spectrum  is 

03  =  c/fc, 


ftw  =  2  [c?  +  c?-^  +  K')] 

+  C1  -  +  K2)]2  -  cfc?[  1  -  p(ic*  +  k2z)  +  4 <pc2(c2  -  cJ)(k2xk2z)}}1/2, 


V 


where 


¥>  = 


A  / 


A  l  +  Ah 


(4  +  *J), 


is  a  kind  of  order  parameter.  The  c*,  q  are  the  usual  transverse  and  longitudinal  sound- 
speeds.  It  is  seen  that  the  branch  corresponding  to  02  has  the  peculiarity  of  disappearing 
at  the  approach  of  the  phase  transition, 


(Ah  <C  1), 

in  the  directions  of  the  magnetic  field  and  of  the  easy  axis.  It  is  the  nature  of  boson 
statistics  that  as  a  quasiphonon  branch  goes  to  zero  in  the  neighborhood  of  the  critical 
field,  the  total  energy  stored  in  the  corresponding  collective  excitations, 

E  = 


increases.  In  other  words,  one  is  not  penalized  by  transferring  bosons  to  this  mode. 
This  means  that  in  the  equilbrium  corresponding  to  the  increased  field,  energy  will  be 
transferred  to  quasimagnons  from  quasiphonons.  This  will  mean  a  lowering  of  the  tem¬ 
perature,  thus  one  would  expect  a  negative  magnetocaloric  effect  in  the  neighborhood 
of  the  phase  transition.  A  similar  effect  wherein  occurs  change  of  polarization  of  sound 
waves  has  been  observed  in  experiments  with  AC  magnetic  pulses.  It  would  appear  that 
a  simplified  model  is  called  for  that  can  be  integrated  over  the  whole  Brilluoin  zone  in 
order  to  compute  the  entropy  but  still  reproduces  the  behavior  of  this  spectrum  along 
the  symmetry  axes.  In  the  simplest  model  of  this  sort  the  spectrum  could  be  written 
in  the  form: 


magnon  frequencies 


T,2k  =  [ ak 2  +  n3][n3cos20  +  ak2], 


phonon  frequencies 


k2{n0  ± 
n2k2, 


—  Til  cos 26~}, 
k 


(20) 

(21) 


where 

cos20  OC  1 1  —  h\  =  j - r—— — |. 

Ha 

The  isentropic  curves  for  this  system  show  a  very  strong  negative  magnetocaloric  effect 
at  a  field  at  a  finite  distance  below  the  critical  field. 


According  to  the  Ginzburg- Landau-Abrikosov  theory,  the  flux  line  lattice  of  vortices 
in  a  type-II  superconductor  has  a  “magnetoelastic”  interaction  built  into  it.  If  we 
assume  a  linear  elastic  description  of  the  flux  line  lattice,  then  there  are  independent 
modes  of  oscillation,  shear,  compression,  and  tilt  waves;  all  of  these  carry  the  magnetic 
interaction  with  them.  For  the  high-temperature  superconductors  with  their  large  values 
of  the  Ginzburg-Landau  parameter,  k,  we  should  expect  a  large  negative  magnetocaloric 
effect.  For  small  values  of  the  flux,  the  elastic  energy  is  proportional  to  the  square  of 
the  flux,  so  we  should  expect  a  positive  magnetocaloric  effect,  as  for  the  ferromagnet, 
only  for  vanishing  fields.  Both  the  flux  penetration  depth,  A,  and  the  coherence  length, 
(,  diverge  as  «  (1  —  5/5c2)~1j/2  as  the  flux  approaches  the  upper  critical  field,  Bc2.  This 
means  a  vanishing  of  the  elastic  exitation  energy  of  the  flux  line  lattice  as  the  upper 
critical  field  is  approached,  leading  us  to  consider  approaching  this  situation  in  a  way 
analogous  to  the  magnon-phonon  interaction.  For  a  uniaxial  flux  line  lattice  where  the 
flux  is  always  perpendicular  to  the  ab-plane,  and  we  ignore  elastic  interactions  with 
pinning  caused  by  material  inhomogeneities,  the  flux-line  lattice  modes  have  the  elastic 
energy  matrix  elements 


<Aj  =  (cn(b,k)  -  <^e,(b))ktk3  +  ^j(c66(6)(fc2  +  k2y)  +  c44{k,b)k2z) 


The  computed  isentropic  curves  for  shear  excitations  show  a  monotone  negative  mag¬ 
netocaloric  effect.  In  the  isentropic  curves  for  excitations  corresponding  to  isotropic 
compression  waves  there  is  an  interesting  change  of  sign  of  the  magnetocaloric  effect 
at  a  low  value  of  the  magnetic  field.  Previously,  such  behavior  had  been  associated 
with  antiferrimagnetic  systems.  We  have  identified  a  possible  test  material,  a  nanocrys¬ 
talline  superparamagnetic  material.  Besides  being  designed  to  match  a  certain  model 
of  magnetoelasticity,  this  material  system  is  an  example  of  a  mesoscopic  system  with 
significant  interactions  in  the  nanometer-micron  range.  The  manipulation  and  process 
control  of  the  spatial  and  temporal  structure  of  matter  over  mesoscopic  distance  scales, 
especially  near  phase  transitions,  has  been  identified  as  important  to  the  long-range 
war-fighting  capabilites  of  the  Air  Force. 
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sions  with  scientists  from  AFRL  Air  Vehicles  Directorate  and  have  generated  research 
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I  have  run  an  Applied  Math/Physics/Modeling  Seminar.  Speakers  that  were  partially 
supported  by  this  grant  were  M.  Gunzburger  of  Iowa  State  (who  spoke  on  nun'ierical 
methods),  S.  Shen  of  NASA/GSFC  (who  spoke  on  data  mining),  and  V.  Rao  of  the 
Swedish  Technical  Institute  (who  spoke  on  non-linear  susceptibility  in  superconductors). 
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C.  Pettit,  and  T.  Knox.  The  PI  has  also  hosted  visitors  from  University  of  Utah,  Iowa 
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zburger  [most  of  expense  paid  for  by  host].  In  October  99  I  travel  to  College  Station, 
Texas,  and  gave  an  invited  talk  at  the  International  Conference  on  Nonlinear  Control 
Theory  [most  of  expense  paid  for  by  host].  In  December  99,  I  traveled  to  Prague,  CR, 
collaborate  with  Z.  Kalva  on  the  magnetocaloric  effect  [most  of  expense  paid  for  by 
host],  I  presented  a  poster  at  the  American  Physical  Society  March  2000  meeting  in 
Minneapolis.  In  October  2000,  I  presented  an  invited  talk  at  the  International  Confer¬ 
ence  on  Control  in  Raleigh,  N.C.;  in  November  2000,  I  gave  a  talk  at  the  APS  Annual 
Meeting  of  the  Division  of  Fluid  Dynamics  in  Washington,  D.C.;  in  March  2001,  I  pre¬ 
sented  a  poster  at  the  Annual  Meeting  of  American  Physical  Society  in  Seattle,  WA.; 
in  April  2001,  I  gave  an  invited  talk  at  the  American  Mathematical  Society  Western 
Regional  Conference  in  Las  Vegas,  NV.  I  presented  a  poster  in  March  02  at  the  Annual 
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at  the  SPIE  High  Power  Laser  Ablation,  in  Taos,  NM.  I  have  attended  workshops  on 
computational  methods  in  in  Detroit,  MI;  Kent,  OH;  Columbus,  OH;  Urbana,  IL.  I 
have  given  presentations  at  AFOSR  Contractors  meetings  in  Stanford  (2000,  2001)  and 
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